Nonlinear left-handed metamaterials 
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We analyze nonlinear properties of microstructured materials with the negative refractive index, 
the so-called left-handed metamaterials. We demonstrate that the hysteresis-type dependence of 
the magnetic permeability on the field intensity allows changing the material properties from left- 
to right-handed and back. Using the finite-difference time-domain simulations, we study the wave 
reflection from a slab of nonlinear left-handed material, and observe the propagation of tempo- 
ral solitons in such materials. We demonstrate also that nonlinear left-handed metamaterials can 
support both TE- and TM-polarized self-trapped localized beams, spatial electromagnetic solitons. 
Such solitons appear as single- and multi-hump beams, being either symmetric or antisymmetric, 
and they can exist due to the hysteresis-type magnetic nonlinearity and the effective domains of 
negative magnetic permeability. 
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I. INTRODUCTION 

Recent theoretical studies 0, 0, and experimental 
results [10 have shown the possibility of creating 
novel types of microstructured materials that demon- 
strate many intriguing properties such as the effect of 
negative refraction. In particular, the composite mate- 
rials created by arrays of wires and split-ring resonators 
were shown to possess a negative real part of the magnetic 
permeability and dielectric permittivity for microwaves. 
These materials are often referred to as left-handed mate- 
rials (LHMs) or materials with negative refraction. Prop- 
erties of the left-handed materials were analyzed theoret- 
ically by Veselago long time ago 0, but such materials 
were demonstrated experimentally only very recently. As 
was shown by Veselago , the left-handed materials pos- 
sess a number of peculiar properties, including negative 
refraction for interface scattering, inverse light pressure, 
reverse Doppler and Vavilov-Cherenkov effects, etc. 

So far, most of the properties of left-handed materi- 
als were studied in the linear regime of wave propagation 
when both magnetic permeability and dielectric permit- 
tivity of the material are assumed to be independent 
on the intensity of the electromagnetic held. However, 
any future effort in creating tunable structures where the 
field intensity changes the transmission properties of the 
composite structure would require the study of nonlinear 
properties of such metamaterials, which may be quite 
unusual. In particular, the recently fabricated metama- 
terials are composed of a mesh of wires and split-ring 
resonators (SRRs). The wires provides negative dielec- 
tric permittivity, while SRRs give negative magnetic per- 
meability. Metamaterials possess left-handed properties 
only in some finite frequency range, which is basically 
determined by the geometry of the structure. The pos- 
sibility to control the effective parameters of the meta- 
material using nonlinearity has recently been suggested 
in Refs. [alii- Importantly, the microscopic electric field 
in the left-handed structure can be much higher than 
the macroscopic electric field carried by the propagat- 



ing wave. This provides a simple physical mechanism for 
enhancing nonlinear effects in left-handed materials. 

In this paper we present a brief overview of nonlinear 
properties of left-handed metamaterials for the example 
of a lattice of the split-ring resonators and wires with 
a nonlinear dielectric. By means of the finite-difference 
time-domain (FDTD) simulations, we study the wave 
scattering on a slab of a nonlinear composite structure. 
We also discuss the structure of electromagnetic solitons 
supported by the nonlinear left-handed materials with 
hysteresis-type nonlinear response. We believe our find- 
ings may stimulate the future experiments in this field, as 
well as the studies of nonlinear effects in photonic crys- 
tals, where the phenomenon of negative refraction is an- 
alyzed now very intensively |10t 1 1 1| . 



II. NONLINEAR RESONANT RESPONSE 

First, we follow the original paper Q and consider 
a two-dimensional composite structure consisting of a 
square lattice of the periodic arrays of conducting wires 
and split-ring resonators (SRR). We assume that the 
unit-cell size d of the structure is much smaller then the 
wavelength of the propagating electromagnetic field and, 
for simplicity, we choose the single-ring geometry of a 
lattice of cylindrical SRRs. The results obtained for this 
case are qualitatively similar to those obtained in the 
more involved cases of double SRRs. This type of mi- 
crostructured materials has recently been suggested and 
built in order to create left-handed metamaterials with 
negative refraction in the microwave region Q . 

The negative real part of the effective dielectric per- 
mittivity of such a composite structure appears due to 
the metallic wires whereas a negative sign of the mag- 
netic permeability becomes possible due to the SRR lat- 
tice. As a result, these materials demonstrate the prop- 
erties of negative refraction in the finite frequency band, 
luq < lo < min(o;p, u>\\ m ), where ojq is the eigenfrequency 
of the SRRs, wii m is the frequency of the longitudinal 
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magnetic plasmon, io v is an effective plasma frequency, 
and u) is the angular frequency of the propagating electro- 
magnetic wave, (£, Ti.) ~ (E, H) exp (iwt). The split-ring 
resonator can be described as an effective LC oscillator 
(see, e.g. Ref. with the capacitance of the SRR gap, 
as well as an effective inductance and resistance. 

Nonlinear response of such a composite structure can 
be characterized by two different contributions. The first 
one is an intensity-dependent part of the effective dielec- 
tric permittivity of the infilling dielectric. For simplicity, 
we may assume that the metallic structure is embedded 
into a nonlinear dielectric with a permittivity that de- 
pends on the intensity of the electric field in a general 
form, en = £_d(|E| 2 ). For results of calculations pre- 
sented below, we take the linear dependence that corre- 
sponds to the Kerr-type nonlinear response. 

The second contribution into the nonlinear properties 
of the composite material comes from the lattice of res- 
onators, since the SRR capacitance (and, therefore, the 
SRR eigenfrequency) depends on the strength of the lo- 
cal electric field in a narrow slot. The intensity of the 
local electric field in the SRR gap depends on the elec- 
tromotive force in the resonator loop, which is induced 
by the magnetic field. Therefore, the effective magnetic 
permeability /j, e ff depends on the macroscopic (average) 
magnetic field H, and this dependence can be found in 
the form 



where 



Flu 2 

Meff(H) = 1 + — TTjr 2 , 



w CWVi( H ) 



(1) 



aJ [nhe D (\E g (H)\ 2 )] 



is the eigenfrequency of oscillations in the presence of the 
external field of a finite amplitude, h is the width of the 
ring, r = c 2 /2%aah, for h < S, and T — c 2 /2iraa8, for 
h > S. It is important to note that Eq. has a simple 
physical interpretation: The resonant frequency of the 
artificial magnetic structure depends on the amplitude 
of the external magnetic field and, in turn, this leads to 
the intensity-dependent function /i e g. 

Figures ^ an d HI summarize different types of nonlin- 
ear magnetic properties of the composite, which are de- 
fined by the dimcnsionlcss frequency of the external field 
Q = lo/loq, for both a focusing [Figs. ^ |2{a,b)] and a 
defocusing [Figs. ^ Gfod)] nonlinearity of the dielectric. 

Due to the high values of the electric field in the slot 
of SRR as well as resonant interaction of the electro- 
magnetic field with the SRR lattice, the characteristic 
magnetic nonlinearity in such structures is much stronger 
then the corresponding electric nonlinearity. Therefore, 
magnetic nonlinearity should dominate in the composite 
metamaterials. More importantly, the nonlinear medium 
can be created by inserting nonlinear elements into the 
slots of SRRs, allowing an easy tuning by an external 
field. 
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FIG. 1: Real part of the effective magnetic permeability vs. 
intensity of the magnetic field: (a) fl > 1, a = 1; (b) fl < 1, 
a = 1, (c) f2 > 1, a = —1; and (d) Q < 1, a = —1. Black 
- the lossless case (7 = 0), grey-the lossy case (7 = 0.05). 
Dashed curves show unstable branches. 
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FIG. 2: Imaginary part of the effective magnetic permeability 
vs. intensity of the magnetic field for 7 = 0.05: (a) fl > 1, 
a = 1; (b) H < 1, a = 1, (c) n > 1, a = -1; and (d) fi < 1, 
a — — 1. Dashed curves show unstable branches. 



The critical fields for switching between LH and RH 
states, shown in the Figs.n can be reduced to a desirable 
value by choosing the frequency close to the resonant 
frequency of SRRs. Even for a relatively large differ- 
ence between the SRR eigenfrequency and the external 
frequency, as we have in Fig. ^b) where SI = 0.8 (i.e. 
lu = 0.8wo) ; the switching amplitude of the magnetic 
field is ~ 0.03_E C . The characteristic values of the focus- 
ing nonlinearity can be estimated for some materials such 
as n-InSb for which E c = 200V /cm [r|. As a result, the 
strength of the critical magnetic field is found as H c \ « 
l.6A/m. Strong defocusing properties for microwave fre- 
quencies are found in Ba^Sri^TiOs (see Ref. 01 an d 
references therein). The critical nonlinear field of a thin 
film of this material is E c = 4 • 10 4 V/ cm, and the corre- 
sponding field of the transition from the LH to RH state 
[see Fig. |T](c)] can be found as H c w 55AA/m. 

The unique possibility of strongly enhanced effective 
nonlinearities in the left-handed metamaterials revealed 
here may lead to an essential revision of the concepts 
based on the linear theory, since the electromagnetic 
waves propagating in such materials always have a finite 
amplitude. At the same time, the engineering of non- 



3 



linear composite materials may open a number of their 
novel applications such as frequency multipliers, beam 
spatial spectrum transformers, switchers, limiters, etc. 



III. FDTD SIMULATIONS OF NONLINEAR 
TRANSMISSION 

In order to verify the specific features of the left- 
handed metamaterials introduced by their nonlinear re- 
sponse, in this section we study the scattering of elec- 
tromagnetic waves from the nonlinear metamaterial dis- 
cussed above. In particular, we perform the FDTD nu- 
merical simulations of the plane wave interaction with a 
slab of LHM of a finite thickness. We use the Maxwell's 
equations in the form 

VxE = --— , 

c at 

1 <9E 47r 

VxB = -l- + - (j)+47rVxM, (2) 
c at c 

where (j) is the current density averaged over the period 
of the unit cell, and M is the magnetization of the meta- 
material. We base our analysis on the microscopic model 
recently discussed in Ref. [lfj , and write the constitutive 
relations in the form 

<rM*$- + 0) = ^E, 



dt 
M 



d? 

17*171 2 t B 

— ™ I R — 



2c 



(3) 



where L w is the inductance of the wire per unit length, 
a is the conductivity of metal in the composite, S is the 
effective cross-section of a wire, S ~ Trr^,, for 8 > r w , 
and S ~ 7T(5(2r — <5), for S < r w , where 5 = cj \J2-kglo is 
the skin-layer thickness, Ir is the current in SRR, n m is 
concentration of SRRs. The current in SRRs is governed 
by the equation 



din 

' dt 



dH 1 



dt 



U - RI R , 



(4) 



where L is inductance of the SRR, R is resistance of 
the SRR wire, U is the voltage on the SRR slit, and H' 
is the acting (microscopic) magnetic field, which differs 
from the average (macroscopic) magnetic field. Voltage 
U at the slit of SRR is coupled to the current / through 
the relation 



(5) 



with 



C(U) = irr 2 e (l + a \U\ 2 /f/ c 2 ) /4ird gi 



where e is the linear part of the permittivity of a dielectric 
material inside the SRR slit, U c is the characteristic non- 
linear voltage, and a = ±1 corresponds to the case of the 
focusing and defocusing nonlinear response, respectively. 
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FIG. 3: Reflection (solid) and transmission (dashed) coeffi- 
cients for a slab of nonlinear metamaterial vs. the incident 
field intensity in a stationary regime, for the case of defo- 
cusing nonlinearity (a = —1). Inset shows real (solid) and 
imaginary (dashed) parts of the magnetic permeability inside 
the slab. 



The microscopic magnetic field H' can be expressed in 
terms of M and B using the Lorenz-Lorentz relation ^(| : 



H' = B 



8tt 



M. 



(6) 



As a result, Eqs. ©to © form a closed system of cou- 
pled equations, and they can be solved numerically using, 
for example, the numerical FDTD method. We also no- 
tice that, by substituting the harmonic fields into these 
equations, we recover the expression for the magnetic 
permeability Q. 

Our goal is to study the temporal dynamics of the wave 
scattering by a finite slab of nonlinear metamaterial. For 
simplicity, we consider a one-dimensional problem that 
describes the interaction of the plane wave incident at 
the normal angle from air on a slab of metamaterial of 
a finite thickness. We consider two types of nonlinear 
effects: (i) nonlinearity-induced suppression of the wave 
transmission when initially transparent left-handed ma- 
terial becomes opaque with the growth of the input am- 
plitude, and (ii) nonlinearity-induced transparency when 
an opaque metamaterial becomes left-handed (and there- 
fore transparent) with the growth of the input amplitude. 
The first case corresponds to the dependence of the ef- 
fective magnetic permeability on the external field shown 
in Figs. ^a,c), when initially negative magnetic perme- 
ability (we consider e < in all frequency range) be- 
comes positive with the growth of the magnetic field in- 
tensity. The second case corresponds to the dependence 
of the magnetic permeability on the external field shown 
in Figs.mb). 

In all numerical simulations, we use linearly growing 
amplitude of the incident field within the first 50 periods, 
that becomes constant afterwards. The slab thickness is 
selected as 1.3A where A is a free-space wavelength. For 
the parameters we have chosen, the metamaterial is left- 
handed in the linear regime for the frequency range from 
fx = 5.787 GHz to f 2 = 6.05 GHz. 

Our simulations show that for the incident wave with 
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Incident field intensity 

FIG. 4: Reflection (solid) and transmission (dashed) coeffi- 
cients for a slab of nonlinear metamaterial vs. the incident 
field intensity in a stationary regime, for the focusing non- 
linearity (q = 1). Inset shows real (solid) and imaginary 
(dashed) parts of the magnetic permeability inside the slab. 
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FIG. 5: (a) Reflected (solid) and incident (dashed) wave 
intensity vs time for small amplitudes of the incident wave (i.e. 
in the linear regime). (b,c) Distribution of the magnetic and 
electric fields at the end of simulation time; the metamaterial 
is shaded. 



the frequency /o = 5.9 GHz (i.e. inside the left- 
handed transmission band) , electromagnetic field reaches 
a steady state independently of the sign of the nonlin- 
earity. Both reflection and transmission coefficients in 
the stationary regime are shown in Figs, l-'il II as func- 
tions of the incident field amplitude, for both defocus- 
ing and focussing nonlinear response of the dielectric in- 
filling the SRR slits. In the linear regime, the effective 
parameters of the metamaterial at the frequency /o are: 
e = —1.33 — O.Oli and fi = —1.27 — 0.3i; this allows ex- 
cellent impedance matching with surrounding air. The 
scattering results in a vanishing reflection coefficient for 
small incident intensities (see Figs. 1314f t . 

Reflection and transmission coefficients are qualita- 
tively different for two different types of infilling nonlin- 
ear dielectric. For the defocusing nonlinearity, the reflec- 
tion coefficient varies from low to high values when the 
incident field exceeds some threshold value (see Fig. 01. 
Such a sharp transition can be explained in terms of the 
hysteresis behavior of the magnetic permeability shown 
in Fig. njc). When the field amplitude in metamate- 
rial becomes higher than the critical amplitude [shown 
by a dashed arrow in Fig. ^c)], magnetic permeability 
changes its sign, and the metamaterial becomes opaque. 
Our FDTD simulations show that for overcritical ampli- 
tudes of the incident field, the opaque region of positive 
magnetic permeability appears inside the slab (see the 
inset in Fig. OJ - The magnetic permeability experiences 
an abrupt change at the boundary between the trans- 
parent and opaque regions. The dependencies shown in 
Fig. |3| are obtained for the case when the incident field 
grows from zero to a steady-state value. However, tak- 
ing different temporal behavior of the incident wave, e.g. 
increasing the amplitude above the threshold value and 
then decreasing it to the steady state, one can get differ- 
ent values of the stationary reflection and transmission 
coefficients, and different distributions of the magnetic 
permeability inside the metamaterial slab. Such prop- 
erties of the nonlinear metamaterial slab are consistent 
with the predicted multi- valued dependence of the mag- 
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FIG. 6: The same as in Fig. 5 but in the regime of overcritical 
amplitude of the incident wave. 



netic permeability on the amplitude of the magnetic field. 

In the case of focussing nonlinearity (see Fig. QJ, the 
dependence of the reflection and transmission coefficients 
on the amplitude of the incident field is smooth. This 
effect originates, firstly, from a gradual detuning from the 
impedance matching condition, and, for higher powers, 
from the appearance of an opaque layer (see the inset in 
Fig-EJ with a positive value of the magnetic permeability 
that is a continuous function of the coordinate inside the 
slab. 

Now we consider another interesting case when ini- 
tially opaque metamaterial becomes transparent with the 
growth of the incident field amplitude. We take the fre- 
quency of the incident field to be /q = 5.67 GHz, so that 
magnetic permeability is positive in the linear regime and 
the metamaterial is opaque. In the case of self-focussing 
nonlinear response (a — 1), it is possible to switch the 
material properties to the regime with negative magnetic 
permeability [see Fig^b) making the material slab left- 
handed and therefore transparent. Moreover, one can 
expect the formation of self-focused localized states in- 
side the composite, the effect which was previously dis- 
cussed for the interaction of the intense electromagnetic 
waves with over-dense plasma |l7j . Figure^a) shows the 
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FIG. 7: The same as in Fig. 5 but in the regime of strongly 
overcritical amplitude of the incident wave. 



temporal evolution of the incident and reflected wave in- 
tensities for small input intensities, this case corresponds 
to the linear regime. The reflection coefficient reaches a 
steady state after approximately 100 periods. The spatial 
distribution of the electric and magnetic fields at the end 
of simulation time is shown in Fig. EIb,c), respectively. 

In a weakly nonlinear overcritical regime (see Fig. [SJl , 
the intensity of the reflected beam decreases approaching 
a steady state. In this case we observe the formation of a 
localized state inside the metamaterial slab and near the 
interface, as can be seen more distinctly in Fig.[|Jc). This 
effect give an additional contribution to the absorption 
of the electromagnetic energy, thus leading to a decay of 
the value of the reflection coefficient. 

In a strongly nonlinear overcritical regime, we observe 
the effect of the dynamical self-modulation of the re- 
flected electromagnetic wave that results from the pe- 
riodic generation of the self-localized states inside the 
metamaterial (see Fig.[7J). Such localized states resemble 
temporal solitons, which transfer the energy away from 
the interface. Figure [7fc) shows an example when two 
localized states enter the metamaterial. These localized 
states appear on the jumps of the magnetic permeabil- 
ity and, as a result, we observe a change of the sign of 
the electric field derivative at the maximum of the soli- 
ton intensity, and subsequent appearance of transparent 
regions in the metamaterial. Unlike all previous cases, 
the field structure in this regime do not reach any steady 
state for high enough intensities of the incident field. 



IV. ELECTROMAGNETIC SPATIAL SOLITONS 

Similar to other nonlinear media [Tsj . nonlinear left- 
handed composite materials can support self-trapped 
electromagnetic waves in the form of spatial solitons. 
Such solitons possess interesting properties because they 
exist in materials with a hysteresis-type (multi-stable) 
nonlinear magnetic response. Below, we describe novel 
and unique types of single- and multi-hump (symmetric, 
antisymmetric, or even asymmetric) backward-wave spa- 
tial electromagnetic solitons supported by the nonlinear 
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FIG. 8: Examples of different types of solitons: (a) fun- 
damental soliton; (b,c) solitons with one domain of negative 
or positive magnetic permeability (shaded), respectively; (d) 
soliton with two different domains (shaded). Insets in (b,c) 
show the magnified regions of the steep change of the mag- 
netic field. 



magnetic permeability. 

Spatially localized TM-polarized waves that are de- 
scribed by one component of the magnetic field and two 
components of the electric field. Monochromatic station- 
ary waves with the magnetic field component H = H y 
propagating along the z-axis and homogeneous in the 
y-direction, [~ exp (iu>t — ikz)], are described by the di- 
mensionless nonlinear Helmholtz equation 



d 2 H 

dx 2 



[e l x{\H\ 2 )-l 2 \H = 0, 



(7) 



where 7 = fc/w is a wavenumber, x = x'lu/c is the di- 
mensionless coordinate, and x' is the dimensional coor- 
dinate. Different types of localized solutions of Eq. Q 
can be analyzed on the phase plane {H,dH/dx) (see, 
e.g., Refs. 18]). First, we find the equilibrium points: 
the point (0, 0) existing for all parameters, and the point 
(0,-ffi), where H\ is found as a solution of the equation 



* 2 (#i) 



Y 2 



CI 2 



1 



il 2 - e c ff) 



(8) 



Below the threshold, i.e. for 7 < 7t r , where j 2 T = 
e[l + Ffl 2 /(1 - ft 2 )}, the only equilibrium state (0,0) 
is a saddle point and, therefore, no finite-amplitude or 
localized waves can exist. Above the threshold value, i.e. 
for 7 > 7 tr , the phase plane has three equilibrium points, 
and a separatrix curve corresponds to a soliton solution. 

In the vicinity of the equilibrium state (0,0), linear 
solutions of Eq. Q describe either exponentially grow- 
ing or exponentially decaying modes. The equilibrium 
state (0,iii) describes a finite-amplitude wave mode of 
the transverse electromagnetic field. In the region of 
multi-stability, the type of the phase trajectories is de- 
fined by the corresponding branch of the multi-valued 
magnetic permeability. Correspondingly, different types 
of the spatial solitons appear when the phase trajecto- 
ries correspond to the different branches of the nonlinear 
magnetic permeability. 
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The fundamental soliton is described by the separa- 
trix trajectory on the plane (H,dH/dx) that starts at 
the point (0,0), goes around the center point (0,-ffi), 
and then returns back; the corresponding soliton profile 
is shown in Fig.[3Ja). More complex solitons are formed 
when the magnetic permeability becomes multi-valued 
and is described by several branches. Then, soliton so- 
lutions are obtained by switching between the separatrix 
trajectories corresponding to different (upper and lower) 
branches of magnetic permeability. Continuity of the 
tangential components of the electric and magnetic fields 
at the boundaries of the domains with different values of 
magnetic permeability implies that both H and dH/dx 
should be continuous. As a result, the transitions be- 
tween different phase trajectories should be continuous. 

Figures|HIb,c) show several examples of the more com- 
plex solitons corresponding to a single jump to the lower 
branch of fi(H ) (dotted) and to the upper branch of ^(H) 
(dashed), respectively. The insets show the magnified 
domains of a steep change of the magnetic field. Both 
the magnetic field and its derivative, proportional to the 
tangential component of the electric field, are continu- 
ous. The shaded areas show the effective domains where 
the value of magnetic permeability changes. Figure |SJd) 
shows an example of more complicated multi-hump soli- 
ton which includes two domains of the effective magnetic 
permeability, one described by the lower branch, and the 
other one - by the upper branch. In a similar way, we 
can find more complicated solitons with different number 
of domains of the effective magnetic permeability. 

We note that some of the phase trajectories have dis- 
continuity of the derivative at H = caused by infinite 
values of the magnetic permeability at the correspond- 
ing branch of /i c ff(H). Such a non- physical effect is an 
artifact of the lossless model of a left-handed nonlinear 
composite considered here for the analysis of the soliton 
solutions. In more realistic models that include losses, 
the region of multi-stability does not extend to the point 
H = 0, and in this limit the magnetic permeability re- 
mains a single-valued function of the magnetic field 

For such a multi- valued nonlinear magnetic response, 
the domains with different values of the magnetic per- 
meability "excited" by the spatial soliton can be viewed 



as effective induced left-handed waveguides which make 
possible the existence of single- and multi-hump soliton 
structures. Due to the existence of such domains, the 
solitons can be not only symmetric, but also antisym- 
metric and even asymmetric. Formally, the size of an ef- 
fective domain can be much smaller than the wavelength 
and, therefore, there exists an applicability limit for the 
obtained results to describe nonlinear waves in realistic 
composite structures. 

When the infilling dielectric of the structure displays 
self-focusing nonlinear response, we have fl < 1, and in 
such system we can find dark solitons, i.e. localized dips 
on the finite-amplitude background wave fl9l | . Similar to 
bright solitons, there exist both fundamental dark soli- 
tons and dark solitons with domains of different values 
of magnetic permeability. For self-defocusing nonlincar- 
ity and fl < 1, magnetic permeability is a single- valued 
function, and such a nonlinear response can support dark 
solitons as well, whereas for self-focusing dielectric, we 
have Q > 1 and no dark solitons can exist. 



V. CONCLUSIONS 

We have discussed novel properties of left-handed 
metamaterials associated with their nonlinear resonant 
response. For the case of harmonic fields, we have calcu- 
lated the effective magnetic permeability of microstruc- 
tured materials consisting of rods and split-ring res- 
onators, and predicted the hysteresis- like dependence of 
the nonlinear magnetic permeability as a function of the 
applied magnetic field. Using the finite-difference time- 
domain numerical simulations, we have studied the tem- 
poral dynamics of the wave reflection from a slab of 
nonlinear metamaterial that is found to be consistent 
with our theory. Finally, we have predicted the exis- 
tence of electromagnetic spatial solitons supported by the 
hysteresis-type nonlinear magnetic permeability of a left- 
handed material. 

The work has been supported by the Australian Re- 
search Council and the Australian National University. 
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